A Magneto-Optical Trap for Polar Molecules 
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We propose a method for laser cooling and trapping a substantial class of polar molecules, and in particular 
titanium (II) oxide (TiO). This method uses pulsed electric fields to nonadiabatically remix the ground-state 
magnetic sublevels of the molecule, allowing us to build a magneto-optical trap (MOT) based on a quasi-cycling 
J' = 7" — 1 transition. Monte-Carlo simulations of this electrostatically remixed MOT (ER-MOT) demonstrate 
the feasibility of cooling TiO to a temperature of 10 and trapping it with a radiation-pumping-limited lifetime 
on the order of 80 ms. 

PACS numbers: 37.10.Pq, 37.10.Mn, 37.10.Vz 



The field of ultracold polar molecules has recently made 
great strides. Coherent optical transfer of magneto-associated 
molecules can now produce ultracold molecular gases in the 
X 1 ^ (v = 0) ground state with densities of 10 12 cm~ 3 and 
translational temperatures of 350 nK fljj] . Incoherent photo- 
association techniques can reach the X L (v = 0) ground state 
at 100 jtiK [2]. With these temperatures and the reasonably 
large electric dipole moments available from heteronuclear 
bialkali molecules (e.g. 0.76 D for X 1 ! (v = 0) KRb H), 
progress towards quantum simulations of condensed matter 
systems [J, |5[] and quantum computation ja M] should be 
rapid. In fields such as ultracold chemistry [8], access to 
molecular species beyond the bialkali family is of great in- 
terest. Arbitrary species can be cooled to the kelvin regime 
through buffer-gas cooling [91 1 lOfl . while Stark deceleration 
| UJ, Ll2[] reaches the tens of millikelvin level for selected light 
molecules. Unfortunately, there is no demonstrated tech- 
nique to further compress and cool the lukewarm molecular 
clouds resulting from the latter two techniques. Even cavity- 
mediated schemes for molecular laser cooling lfl3i[l^[l5l[l6ll . 
while in the abstract highly attractive methods for cooling 
a broad, chemically interesting set of molecules, have so 
far been unable to cool these lukewarm samples, due to the 
schemes' low scattering rates [15], small cavity mode vol- 
umes 1 1611. and requirement of multiparticle collective effects 

Direct, free-space laser cooling and trapping would be the 
ideal method for producing ultracold molecules, just as it is 
for atoms. Unfortunately, atoms are in general much easier to 
laser cool than molecules, due to the latter's glaring lack of cy- 
cling transitions. Laser cooling generally requires electronic 
transitions, as vibrational and rotational transitions have im- 
practically long excited state lifetimes unless a cavity is used 
11611 . Unfortunately, these "electronic" transitions are never 
purely electronic. Rather, they are rovibronic, and decay into 
various rotational, vibrational, or hyperfine excited states, as 
well as the original ground state 11711 . 

The branching ratios of these rovibronic decays, however, 
are governed by the molecular structure and the dipole selec- 
tion rules. This implies that a clever choice of molecule can 
greatly reduce the number of possible decays. Decays into 
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FIG. 1: (color online, not to scale), (a) The electronic level struc- 
ture of TiO and the transitions of interest for laser cooling. The 
X 3 A ground state is split by the spin-orbit interaction into the three 
X 3 Ai_3 sublevels, of which the X 3 A; level is the lowest. Each sub- 
level contains a vibrational ladder, while each vibrational level con- 
tains a ladder of rotationally excited states (not shown). 48 Ti 16 
has zero nuclear spin and thus there is no hyperfine structure. The 
ground-state A doublet (not shown) is much less than the natural 
linewidth of the E 3 n <— X 3 A transition. The solid arrow denotes the 
v — <— v" = P(l)-branch cooling laser, and the dashed arrows 
denote the v = <— v = 1 and v = <— v = 2 P(l)-branch repump 
lasers. The squiggly lines depict the dipole-allowed decays, with the 
associated Franck-Condon factor q 1 19] next to each decay, (b) The 
rotational and A— doublet structure of the E 3 TI() electronic excited 
state. The states are interleaved, as the rotational splitting is smaller 
than the A— doublet splitting; a and b denote the parity states. Both 
the cooling and repump lasers address the J' = 0, a state. 



excited hyperfine states are impossible in molecules with zero 
nuclear spin, as these molecules have no hyperfine structure. 
The Franck-Condon ratios for decay back to the ground vibra- 
tional level can be quite good (99+% for selected molecules 
1 1811). However, the only constraint on decays to rotationally 
excited levels is that all decays satisfy the total angular mo- 
mentum selection rule AJ = 0,±1 [11]. Thus, even without 
hyperfine structure, it may take up to three lasers per vibra- 
tional level to repump the three possible rotational decays. 
However, if the ground state angular momentum /" is greater 
than the excited state angular momentum /', two of these three 
decays are forbidden. In this case, the molecule must fol- 
low the angular momentum cycle J" —>J' = J" — l—> J", and 
so only one laser is required per relevant vibrational level — 
making laser cooling of these molecules truly practical. 
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Thus, by combining these various transition closure crite- 
ria, we can identify a class of molecules that are exceptionally 
good candidates for laser cooling: they have no net nuclear 
spin, good Franck-Condon overlaps, and their ground or low- 
est metastable state has a higher angular momentum than the 
first accessible electronically excited state. For non-singlet 
molecules, the excited electronic level must also not be a E 
state, as the lack of spin-orbit splitting in E states means that 
the excited state can decay across the spin-orbit ladder. 

We have identified a number of molecules that satisfy all 
of the above requirements. TiO and TiS are both satisfactory 
in their absolute ground states. Metastable FeC, ZrO, HfO, 
ThO, SeO, and the like are promising lf20ll . We expect that 
some other, as-yet uncharacterized metal oxides, sulfides, and 
carbides should also have the necessary electronic structure. If 
one is willing to accept some hyperfine structure as the price 
of chemical diversity, some metal hydrides and metal halides 
may provide additional suitable candidates. 

Of these candidates, we chose to focus on TiO, due to its 
viability in its absolute ground state and the breadth of sp ec- 
troscopy and theory available in the literature 



TABLE I: The wavelengths, Franck-Condon factors, and saturation 
intensities of the cooling and repump transitions of TiO. 
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261 1271 128H . A simplified level structure of TiO is 
shown in Fig.[T] The lowest ground state of TiO is the X 3 Ai, 
with spin-orbit constant A' x a ' = 50.61 cm" 1 and rotational 
constant B (x A > =0.534 cm B28I1 . The lowest excited state is 
theE 3 n withA( E3n ) = 86.82cm- 1 mdB^ =0.515cm- 1 
|2J. (While the d [ E+ and a'A level are energetically below 
the E 3 n level, the inter-system branching ratio is expected to 
be very small.) As Fig. Q] shows, the Franck-Condon fac- 
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21ll for the E 3 F1 — X 3 A band are quite favorable, 
yielding a population leak of 3 x 10 -4 scatter -1 (or a mean 
of ~ 3300 scatters before going into a dark state) with two 
repump lasers. The laser wavelengths, saturation intensities 
Csat)> an d Franck-Condon factors for the cooling and repump 
lines are summarized in TableU Note that these transitions are 
all accessible with diode lasers. The saturation intensities are 
extremely low, as the natural linewidth y of the E 3 n — X 3 A 
transition is on the order of 2% x 32-40 kHz JH Q. This 
is about 5 times weaker than the intercombination line used 
to build a Yb magneto-optical trap (MOT) in 1 2911 . However, 
with the use of a cryogenic buffer-gas-cooled TiO source (sim- 
ilar to [30]), the scattering rate is still large enough to work 
with. 

The prospect of building a TiO MOT is tantalizing, given 
this quasi-closed transition. Traditional MOTs work using a 
/' = J" + 1 transition and a magnetic field to break the degen- 
eracy between the excited-state magnetic sublevels. The MOT 
beams are polarized so that the local orientation and strength 
of the quadrupole magnetic field causes the atom to pref- 
erentially scatter from the laser beam providing a position- 
dependent restoring force and a velocity-dependent damping 
force. The fact that J' > J" means that the atom can always 
scatter from the correct beam, as shown in Fig. |2ja). The 
standard MOT will therefore not work for molecules using 
the aforementioned J' = J" — 1 transition [Fig. |2jb)]. While 
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a magnetic field can break the degeneracy of the ground- 
state magnetic sublevels and thus provide beam selectivity, the 
\nijn\ = J" stretched states are effectively dark states, as they 
can only interact with one of the laser beams, not both lf3ltl . 

What is needed, then, is a way to continually remix the 
ground-state sublevels so that all the molecules spend some 
fraction of their time in bright states. Fortunately, polar 
molecules provide just the handle needed to accomplish this: 
the effective magnetic (B) and electric (E) moments of a polar 
molecule depend in different ways on mj. Thus, applying a 
sudden (i.e. nonadiabatic) electric field orthogonal (or at least 
nonparallel) to the local magnetic field reprojects the total an- 
gular momentum against a new axis, randomizing mj (and 
the A— doublet state) by coupling the two A— doublet mani- 
folds together. At high remix rates and high laser saturations, 
the molecules' time is equally divided across the 2 (2/" + 1) 
ground and 27' + 1 excited states (the factor of 2 in the ground 
state is due to the electrostatic mixing of the A doublet), but 
they can only decay while they are in an excited state. Thus, 
while the molecules are effectively always bright, the maxi- 
mum photon scattering rate is only 2{2j"+i)+(2j'+i) ^' Such 
remixing of the ground-state magnetic sublevels allows the 
building of a new kind of trap, the electrostatically remixed 
magneto-optical trap (ER-MOT). The ER-MOT operation is 
shown in Fig. |2fb). Note that, as the local direction of the 
quadrupole B-field spans all of An steradians over the MOT 
volume, a single E-field pulse will be parallel to the local B- 
field in some region and therefore ineffective at remixing the 
mj's there. This hole can, however, easily be closed by ap- 
plying a second E-field pulse, nonparallel to the first. A basic 
ER-MOT design is shown in Fig.|2c). 

To build an ER-MOT with TiO, there is a minor technical 
complication. To leading order, the molecular magnetic mo- 
ment can be written as fx = jigm j (g^A + g^E) Since 

gL ~ 1 and gs ~ 2, the magnetic moment of the X 3 Ai (A = 2 
and E = — 1) state is small, likely on the order of afls, or the 
fine-structure constant times the Bohr magneton. In contrast, 
while Q. = in the E 3 IIo state, the large A— doublet splitting 
indicates strong mixing with higher electronic excited states, 
and so by analogy with the B 3 IIo optical Zeeman measure- 
ments of I32I1 . we estimate the magnetic moment to be ~ 10 
times that of the X 3 A]. This, combined with the narrowness 
of the E 3 II <-> X 3 A transition, implies that the dynamics of a 
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FIG. 2: (color online), (a) The level structure of a traditional MOT. 
The local magnetic field strength and orientation combined with 
the Doppler shift enhance the scattering from the laser beam that 
provides the damping and restoring forces and suppresses scatter- 
ing from the counter-propagating beam. Since J' > J" , the ground 
state(s) are always able to scatter from every beam, (b) The level 
structure of the ER-MOT. The local magnetic field still governs 
which laser is preferentially scattered, but angular momentum con- 
servation forbids some ground states from interacting with the pre- 
ferred beam. To overcome this, the ground-state magnetic sublevel 
populations are remixed by pulsed electric fields, as represented by 
the dashed lines, (c) A sample ER-MOT design. A pair of elec- 
tromagnet coils are aligned in anti-Helmholtz fashion to produce a 
quadrupole field. Six beams of the cooling laser are converged on 
the center with their polarizations oriented as usual for a MOT, but a 
set of four open-mesh grids are added. The grids are pulsed in pairs 
(e. g., first the X-axis pair and then the Y-axis pair) to produce the 
dipole electric fields needed to remix the magnetic sublevels. The 
center is also illuminated by the repump lasers (not shown). 



TiO ER-MOT will have more in common with narrow-line 
alkaline-earth MOTs [33] than normal alkali metal MOTs. 
Given these predicted magnetic moments, the magnetic gradi- 
ent in a TiO ER-MOT must be < 1 00 G/ cm. This gradient can 
be easily achieved with water-cooled electromagnets Il2ll or 
rare-earth permanent magnets [34]. In contrast, the large (« 3 
Debye [27]) electric dipole moment of TiO and its extremely 
small ground-state A— doublet spacing B26I1 mean that electric 
fields of only IV/ cm will give Stark shifts of about 50y — 
far more than the Zeeman shift within the ER-MOT and thus 
sufficient to reproject mj. These small fields can easily be 
switched with rise times on the order of 10 ns (a frequency of 
2800/ and 56 times the Larmor frequency due to the electric 
field), and thus nonadiabaticity is assured. 

A final concern regarding the viability of the TiO ER-MOT 
is that either the electrostatic or magnetic fields might some- 
how cause population loss by mixing in rotationally excited 
J' > states, which could then decay to J" > 1 states and 
be lost. Fortunately, this loss is inhibited by the ~ 1 cm 
rotational splitting [Fig. [Tib)]. Neither the Zeeman or Stark 
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FIG. 3: (color online). Number (upper, black) and temperature 
(lower, red) time-of-flight plots for the loading of a molecular packet 
into a simulated TiO ER-MOT. The initial spike on the number plot is 
the molecular packet flying through the ER-MOT volume; the broad 
hump is the actual captured molecules. The decay of the molecule 
number is due to radiation pumping of the captured population into 
excited v" > 3 levels, and yields an ER-MOT lifetime of 80 ± 5 ms. 
Error bars represent statistics over multiple simulation runs. 



shifts within the ER-MOT are anticipated to be larger than 
~ 100 MHz, and so the perturbative probability to leave the 

w 10~ 5 . This is much 
1 loss rate from decays to 
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desired J' = state is 
smaller than the 3 x 10~ 4 scatter" 
v" > 3 and so is of no importance. 

To verify the feasibility of building an ER-MOT with TiO, 
we performed a set of 3-D semiclassical Monte-Carlo sim- 
ulations. We conservatively assumed a natural linewidth of 
y = 2n x 32 kHz, a magnetic dipole moment of a/Xg, and an 
electric dipole moment of 3 D. We used a 1 je 1 laser waist 
diameter of 6 cm. Our code treated photon scattering and 
molecule kinematics semiclassically and approximated the 
electrostatic remixing as a sudden reprojection of the diag- 
onalized Zeeman Hamiltonian wavefunction against a new 
Stark + Zeeman Hamiltonian. In addition, the simulations 
used a set of 60 additional red-detuned frequency components 
(4.1y spacing, 7.8 MHz total bandwidth) within the cooling 
beams to increase the capture velocity of the ER-MOT, simi- 
lar to the approach used in the Yb MOT of B29I1 . For ease of 
simulation, we used a discrete molecular packet rather than a 
continuous source. The packet was a 3 mm sphere initially 
centered at —4.2 cm on the X-axis in the coordinate system 
of Fig. I2c). It had a flat velocity distribution of 8 ± 3.5 m/s 
centered in orientation around x, with an opening half-angle 
of 0.3 rad. The magnetic field gradient was 51 G/cm and the 
electric field was pulsed to 4V/cm for 100 ns at a rate of 
50 kHz. The cooling laser was detuned by 3.5/ to the red, 
with a saturation parameter of s zoo \ = 10.5 per frequency com- 
ponent. The repump lasers had saturations 5 le pump =6.1. 

These simulation parameters yield the loading and temper- 
ature curves shown in Fig. |3ja) and Ob), respectively. The 
final temperature was approximately 10 /iK. The ER-MOT 
lifetime was limited to about 80 ms by radiation pumping into 
vibrationally-excited dark states. The estimated capture ve- 
locity was 5.7 m/s, which when taking the rotational distri- 
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FIG. 4: Fractional capture vs electrostatic remix rate after 160 ms of 
simulation time for the TiO ER-MOT. Error bars represent statistics 
over multiple simulation runs. The curve is only a guide to the eye. 

bution into account allows the capture of about 0.02% of a 
4.2 K thermal distribution. If one assumes a 4.2 K source 
flux of 10 10 s _1 [30], the calculated capture velocity, lifetime, 
and ER-MOT radius predict an ER-MOT number of 10 5 and 
a density of 10 9 cm~ 3 . 

In addition, we studied the importance of the electrostatic 
remixing to the ER-MOT operation. Figure [4] plots the 
molecule number after 160 ms against the electrostatic remix 
frequency r rem i x . The plot clearly shows the importance of 
the electrostatic remixing. At r rem i x <C no molecules are 
captured, since the molecules are optically pumped out of the 
bright state much faster than they are remixed back into it. 
For r rem ix < tne ca P ture efficiency rises with increasing 
r rem i x , and then the efficiency saturates around r rem j x = 
as the molecules become evenly divided among the various 
ground sublevels. As additional validation checks on the sim- 
ulation code, we verified that turning off the repump lasers 
does indeed inhibit the formation of an ER-MOT by pumping 
the entire population into vibrationally excited states. We also 
verified that we could reproduce the experimental Yb MOT 
oflH by modifying the code to simulate a J 1 = J" + 1 setup 
with the correct atomic parameters. 

In summary, we have shown that molecules whose lowest 
ground- or metastable level-intersecting electronic transition 
has J ! = J" — 1 constitute good candidates for direct laser 
cooling. We have found several molecules that satisfy this 
requirement and have no hyperfine structure. We have pro- 
posed a method to use the electric dipole moment of these 
molecules to remove dark states in the ground and thus build 



an electrostatically remixed magneto-optical trap. We have 
validated these ideas through Monte-Carlo simulation for a 
specific molecule, TiO, and verified the necessity and efficacy 
of the electrostatic remixing. 
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